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ABSTRACT 

We solve for the evolution of the vertical extent of the convective region of a neutron star atmosphere during a 
Type I X-ray burst. The convective region is well-mixed with ashes of nuclear burning due to the short turbulent 
mixing time scale and its extent determines the rise time of the burst light curve. Using a full nuclear reaction 
network, we show that the maximum vertical extent of the convective region during photospheric radius ex- 
pansion (RE) bursts can be sufficiently great that: (1) some ashes of burning are ejected by the radiation driven 
wind during the RE phase and, (2) some ashes of burning are exposed at the neutron star surface following the 
RE phase. Depending on the ignition conditions, ashes with mass number in the range A 30-60 are mixed in 
with the ejected material. As the ejected material cools during the RE phase some of the ejected heavy-element 
ashes cease to be fully ionized. In addition, those ashes that remain bound to the neutron star will temporarily 
reside in the photosphere after it has settled back down to the neutron star surface. Some of these surface ashes 
are of high enough proton number Z that they are not fully ionized. We calculate the expected column density 
of ejected and surface ashes in hydrogen-like states and determine the equivalent widths of the resulting pho- 
toionization edges from both the wind and neutron star surface. We find that these can exceed 100 eV and are 
potentially detectable. A detection would probe the nuclear burning processes and might enable a measurement 
of the gravitational redshift of the neutron star In addition, we find that in bursts with pure helium burning 
layers, protons from {a, p) reactions cause a rapid onset of the '^C(p,7)'^N(Q:,p)'^0 reaction sequence. The 
sequence bypasses the relatively slow '^C(a, 7)"'0 reaction and leads to a sudden surge in energy production 
that is directly observable as a rapid (s^ ms) increase in flux during burst rise. 

Subject headings: accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances — stars: neu- 
tron — X-rays: bursts 



L INTRODUCTION 

Type I X-ray bursts are produced by the unstable nuclear 
burning of freshly accreted hydrogen- and/or helium-rich ma- 
terial on the surface of a n eutron star (NS) in a low-mass X-ray 
binary (for reviews, s ee iLewin. van Par adiis. & Taaml lT995t 
Hildsten 1998; Strohmayer & Birdsten 2005). The burst en- 
ergies (lO-'^- 10"*" ergs), durations (~ 10- 100 s), and recur- 
rence times (hours to days), depend strongly on the compo- 
sition of the accreted matter and on the accretion rate, M, 
which can range from 10~" to 10~'^Mq yr"'. The burst prop- 
erties are also sensitive to the composition of the ashes of 
burning from preceding bursts ('Taam"1980") as underscored 
by recent burst simulations that impleme nt large nuclear re- 
action netw orks for t he energy generation JSchatz et al.l200lt 
iBrown et al.>,2002; .Wooslev et al..2004) . 

The sensitivity of the nuclear energy generation rate to tem- 
perature and density concentrates the burning during a burst 
to a thin layer at the base of the accreted atmosphere (e.g ., 
iFuiimoto. Hanawa. & Miv aii 1981: iFushiki & LambllT98% . 
Since the thermal time scale there is ~ 1 - 10 s while the dy- 
namical time scale is ~ 10"'' s, the temperature gradient near 
the burning layer is nearly adiabatic resulting in a region of 
highly efficient convection. The short mixing time scale en- 
sures that the ashes of burning are well-mixed throughout the 
convective region. 

The vertical extent of the convective region evolves dur- 
ing the burst, as demonstrated in time-dependent numerical 
simulations (e.g.. Joss 1978; Taam 1980; Avasli & Joss 1982; 
IW)oslev & Weaverill984IWooslev et al.l2004l) . In this paper. 



we carry out a thorough survey of the dependence of the con- 
vection region's extent on M, the composition of the accreted 
material, and the pre-burst thermal state of the atmosphere. 
We also show how the evolution of the convective extent in- 
fluences the observed burst rise times. 

We demonstrate that for photospheric radius expan- 
sion (RE) bursts, in which the super-Eddington luminos- 
ity drives a radiation-driven wind, the convective region 
extends out to sufficiently low pressures that ashes can 
be ejected by the wind. Depending on the burst param- 
eters, the wind can eject as much as ^1% of the ac- 
creted mass ( i Pa^ zynski ^ Pr9SZYn^kil Il986t Lloss & Melial 
ll987HNobiU. TuroUa. & Lapidusll994 . corresponding to flie 
ratio of nuclear energy release to gravitational binding en- 
ergy ~ £nuc/£^srav — 5 McV nuclcon"' /200 MeV nucleon"'. 
Sugimoto. Ebisuzaki. & Hanawa ( 1984) pointed out that such 
mass ejection could expose the underlying helium-rich mate- 
rial and result in an Eddington-limited peak flux for helium 
rather than hydrogen. We show that the ejected ashes may be 
directly detectable with spectroscopy during the RE phase and 
afterwards, when the photosphere, laced with heavy element 
ashes, settles to the NS surface. 

The convective region evolves during a burst in the follow- 
ing sequence of events. As the base temperature rises and the 
nuclear energy generation increases during the early stages 
of a burst, the entropy in the convective region increases. 
Initially, only a negligible amount of thermal energy is lost 
to radiation diffusing into the overlying radiative region and 
underlying crust. Since the time scale for radiative diffu- 
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sion across the convective-radiative interface is longer than 
the burning time scale during these early stages, the con- 
vective region extends vertically outward to lower pressures 
(jpuiimoto, Hanawa, & Mivaii 1981). Eventually, the burning 
rate becomes sufficiently slow that the nuclear energy flux is 
carried most efficiently by radiation rather than by convection. 
At that point the convective region recedes back to higher 
pressures. 

Regardless of whether the convective region is expanding 
outwards or receding downwards, its extent at a given time is 
set by the radial location in the atmosphere where the constant 
entropy of the convective region equals the radially increasing 
entropy of th e overlying radiative region. Based on this ar- 
gument, IJossI ^911) showed that the convective region never 
acquires a high enough entropy to reach the photosphere lo- 
cated at column depth ^ 1 g cm"^. Nonetheless, the convec- 
tive region can reach pressures 5, 1% of that at the base of the 
accreted layer ( Jossl ll978HTaam lll98H: 'Avasli & Jos slll982t 
[Hanawa & Sugimotol [1982t |H anawa &. Fuiimoto 198'^1 In 
their simulations, Wooslev et al. (2004) found that in a burst 
with a pure He burning layer, the peak flux exceeded the Ed- 
dington limit and the convective region extended beyond their 
numerical surface (i.e., the resolution limit of their grid lo- 
cated at a pressure k, 0.3% that at the base). 

In this paper, we solve the time-dependent entropy equa- 
tion that describes the evolving thermal structure of the at- 
mosphere and growth of the convective region. We calculate 
the minimum pressure reached by the convective region (or, 
equivalently, the minimum column depth yc.min) for a range of 
burst parameters. We show that yc,mm is sensitive to the burst 
ignition conditions and that, in general, the larger the burst 
peak flux and the smaller the entropy of the pre-burst atmo- 
sphere, the smaller jcmin- Thus, yc^mm decreases with decreas- 
ing M. We find that in systems where t he accreted m aterial 
is helium-rich, such as 4U 1820-30 (see Cummina 2003 and 
references therein), or in systems accreting solar abundances 
at low M (< 10"^Mo yr"'), jcmin ^ Jwind during RE bursts, 
where ywind is the column depth above which mass is ejected 
by the radiation-driven wind. As a result, the wind ejects some 
of the nuclear burning ashes. Furthermore, some of the ashes 
that remain bound to the NS are exposed at the photosphere. 

In § IJlwe describe our analytic prescription for calculating 
the evolution of the thermal structure of the radiative and con- 
vective regions during an X-ray burst. We explain why the 
evolution during the early, convective stages is sensitive to M 
and accreting composition. In §|3]we describe the full nuclear 
reaction network we use to calculate the nucleosynthesis and 
nuclear energy generation. In §0]we examine how Jcmm de- 
pends on the burst ignition conditions and explore those con- 
ditions most conducive to ash ejection and exposure. 

Since the rise time of the burst light curve is determined by 
the thermal diffusion time at the top of the convective zone, 
we also evaluate the rise time dependence on burst parameters 
such as M and accreting composition. We show these results 
in§El 

In § |5] we show the ash composition profile and discuss 
the observational consequences of heavy-element ash ejec- 
tion and surface exposure. Depending on the ignition con- 
ditions, nuclei as heavy as A 60 are ejected by the wind. 
We calculate the column density of ejected and surface ashes 
in hydrogen-like states and discuss the prospects for detect- 
ing the resulting photoionization edge features during high 
spectral resolution observations of RE bursts. Such features 



probe the nuclear burning and may constrain the NS equa- 
tio n of state. If bursts ig nite in the ashes of previous bursts, 
as 'Wooslev et alJ J200 4V suggest, even heavier elements are 
ejected. These may include some light p-nuclei (Schatz et al] 
1998, 2001), whose origins are not understood. We conclude 
in §|6lwith a summary of our work and mention the possibihty 
of ash ejection during superbursts. 

2. EVOLUTION OF THE ATMOSPHERE DURING A BURST 

In this section we consider the temporal evolution of the 
atmosphere's thermal structure during a burst. We start in 
§Owith a general description of the atmosphere's structure 
and explain the boundary condition at the convective-radiative 
interface. The evolution of the location of the convective- 
radiative interface is described in § 12.21 

2.1. Thermal structure of the atmosphere 

The NS atmosphere maintains hydrostatic equilibrium 
throughout the burst so that pressure varies with height as 
dP/dr = -pg, where p is the density and g is the surface grav- 
ity. We assume a NS mass M = 1 AMq and radius R=10 km, 
giving g = (l+z)GM/R^ = 2.43 x 10'"* cm s"^, where the grav- 
itational redshift z = (1 - IGM/Rc^T^I^ -1=0.31. Since the 
atmosphere is thin compared with the NS radius, g is effec- 
tively constant throughout the accreted layer. Hydrostatic bal- 
ance therefore yields P = gy, where the column depth y, de- 
fined by dy = -pdr, is a convenient parameterization of the 
vertical spatial coordinate. We determine the extent of the 
convective region over a broad range of burst parameters and 
thus consider only one-dimensional models in our calcula- 
tions. We do not account for the affect of a spreading burning 
front during burst rise nor the influence of rotation on the con- 
vective structure, though these effects may be important (see 
Spitkovskv. Levin. & Ushomirskv 2002). 

The entire atmosphere is radiative before helium burning 
at the base triggers the burst. The thermal profile is then de- 
scribed by the diffusion equation 



where F is the outward heat flux and the opacity k has contri- 
butions from electron scattering and free-free absorpti on and 
is calculated using the approximation given by S chatz et alJ 
( 1999). The pre-burst flux Fq = Fn+Fcmst, where Fh is the 
flux from stable hydrogen burning via the hot CNO cycle and 
^cnjst is the flux from heat released by ele ctron captures and 
pycno nuclear reaction s deep in the crust jBrown & BildsteJ 
119981 >Brown 2Q0(k 12004 1). Following the burst igni- 
tion calculations of Gumming (2003, hereafter C03), Fh = 
e//min[yH,3'i] where e// = 5.8 x lO'-' ergs g"' s"'(Z/0.01) is 
the hot CNO energy production rate for a CNO mass frac- 
tion Z, yn is the column depth of the layer that is burning 
hydrogen, and yt, is the column depth at the base. For typ- 
ical ignition conditions y;, ~ 3 x 10^ g cm"^. For a given 
local accretion rate m (in units of g cm"^ s"') and accreted 
hydrogen mass fraction Xy, the hydrogen burning depth is 
yn = 6.8 x lO'^ g cm-2(m/0.1mEdd)(0.01/Z)(Xo/0.71). Here 
mEdd = 2mpc/(l +XQ)Ran = 8.8 x lO"* g cm'^ s-i(1.71/[l + 
Xq]) is the local Eddington accretion rate where mp is the 
proton mass, c the speed of light, and crxh the Thomson 
scattering cross section. For m smaller than the critical ac- 
cretion rate man ~ 0.04;7iEdd, yH < yb and there is enough 
time to burn all the hydrogen at the base before the he- 
lium burning becomes unstable. The burst then ignites in a 
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pure helium layer. As in C03, we assume Fcmst = wiQcmst = 
10^' ergs cm"^ s'^nUQo.i where nu = m/\(f g cm"^ s"' and 
Gciust = O.lQo.i MeV nucleon"' ~ Qo.i 10'^ ergs g"' is the en- 
ergy per nucleon released in the crust from pycnonuclear and 
electron capture reactions that escapes from the surface. The 
pre-burst flux F21 = Fq/IO^' ergs cm"^ s"' is thus, 

F21 =m4Qo.i+mm[37Xo{l+Xo)m,6Zo,oiy&] , (2) 
where Zo.oi =Z/0.01 andyg =yi,/lQ^ g cm"^. 

The thermal evolution of the NS atmosphere during a burst 
is described by the entropy equation 



ds dF 
T— = — + e, 
dt dy 



(3) 



where e is the energy release rate from nuclear burning. Dur- 
ing the burst the entropy grows with time due to nuclear burn- 
ing and we neglect the advective accretion flow. Therefore, 
Tds = CpdT where Cp is the specific heat at constant pres- 
sure. Integrating equation Q over column depth then gives 



TABLE 1 
Burst ignition models 



Model 


ys 


^0 


ZcNO 


Gcriist 


network 


He[;ii] 


3.0 


0.0 


0.01 


0.1 


full 


HeO.la 


3.0 


0.0 


0.01 


0.1 


a-only 


HeO.lb 


5.0 


0.0 


0.01 


0.1 


full 


HeO.lQ 


3.0 


0.0 


0.01 


0.2 


full 


HHe[m] 


3.0 


0.71 


0.01 


0.1 


full 


HHeO.OlX 


3.0 


0.1 


0.01 


0.1 


full 


HHeO.OlXZ 


3.0 


0.1 


0.0001 


0.1 


full 



Note. — Col. (1): Model name, where [m] denote.s 
the accretion rate in units of mndd ■ We consider models 
spanning the range O.Objigdd to 0.2f?iE(id- Col. (2): Igni- 
tion column depth vg = yb/i-0^ gcm"^. Col. (3): Accreted 
hydrogen fraction Xq. Col. (4): CNO mass fraction Zcno- 
Col. (5): Qcru.st in units of MeV nucleon"' . Col. (6): Re- 
action network used. 
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We assume the atmosphere is composed of two regions: a 
completely convective region between yc <y <yb7 and a com- 
pletely radiative region for y < y^. During the burst rise, y/, is 
constant while y^ evolves from an initial value yc = to a min- 
imal value yc = yc,min and finally back to yc = yt- Demarcating 
the atmosphere in this way is reasonable given that the con- 
vective eddies are highly subsonic over most of the convective 
zone, i.e., near the base Vconv — (F /pY^^ ~ 10^ cm s"' ^ Cs ~ 
igyt/py^'^ 2 X 10^ cm s"^ Near the top of the convective 
zone, convection becomes inefficient and Vconv ~ Cj. 

The thermal profile in the radiative region satisfies equa- 
tion Q. Since the radiative region is composed primarily 
of freshly accreted hydrogen and/or helium, the main opac- 
ity is Thomson scattering k ~ ~ (TT:ii{l+X)/2mp (there 
are corrections to Kgs due to relativistic electrons and degen- 
eracy). The opacity varies only slightly with column depth so 
that over much of the radiative region dlnT /d Iny ~ 1 /4. For 
mixed hydrogen/helium accretion at m < merit, a pure helium 
layer develops in the region < y < Ji- In this region Fh = Q 
and since Fcmst is small at low m, the pre-burst profile there is 
nearly isothermal. 

Although the entropy in the convective region grows with 
time, at a given instant it is nearly spatially constant. This, 
in addition to the subsonic motion of the convective ed- 
dies, suggests that the thermal profile in the convective re- 
gion very nearly follows an adiabat so that (dlnT /dlny)com = 
{dlnT/dlnyU = n(y\ i.e., T(yc < y < yt) = Tbiy/ybT^K 
where 7}, = 7i,(f) is the temperature at the base and the adi- 
abatic index n{y) varies with column depth. We define the 
column depth of the convective-radiative interface yc as the 
location where the density of the radiative solution just ex- 
ceeds that of the convective solution (i.e, neutral buoyancy 
criterion). 

For the eq uation of state we use the interpolation formulae 
of lPacz vnski ( 1983) to account for partially degenerate elec- 
trons. Using his notation, the specific heat and adiabatic index 
are given by 
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where P,, P,-, and Pe = (Pgnd +^ed)'^^ are the pressure due 
to ions, radiation, and electrons, respectively, fend and fed 
are an approximation to the degenerate and non-degenerate 
components of the electron pressure, P = P, + Pc + P,-, f = 
dlnPf-ii/dlnp, and 
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At burst onset the pressure is nearly that of an ideal gas and 
n ~ 2/5 while at late times radiation pressure contributes sig- 
nificantly, which in the limit P = Pr gives n= 1/4. 

2.2. Temporal evolution of the thermal structure 

The evolution of the convective extent yc depends on the 
rate at which the base temperature rises, dTh/dt, and the rate 
at which the thermal energy of the overlying radiative region 
increases. The rate of temperature change in the convective 
region is dT /dt = (y/yb)"[dTh/dt + Thln(y/yb)dn/dt]. The 
second term is negligible compared with the first term so that 
by equation Q, 



dt 



r 



edy + FH + ^Lmst - ^ios.s(3'c) 



^Cpiy/ybT^Uy 

'' yc 



(9) 



where FjossCye) is the radiative flux escaping from the convec- 
tive region into the overlying radiative region. Physically, the 
rate at which 7}, changes is determined by the competition be- 
tween the net energy input into the convective region (i.e., the 
energy generated by nuclear burning and crustal heating mi- 
nus the energy lost to radiation) and the energy expended in 
heating up the growing convective region. We describe our 
full nuclear reaction network in §|3] 

We use mixing length theory to estimate Fioss(3't), which 
gives the fraction of the total flux transporte d by convection 
at column depth y (iHansen & KawaleJ 19941) . 
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where V^d = "(y), V,ad = IfnFy j'XacT^ , and Nu = Nu(y) is 
the Nusselt number describing the efficiency of convection. 
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Base Temperature (10^ K) 

Fig. 1. — Thermal timescale at the convective-radiative interface ?thO'c) 
{solid line) and the growth timescale of the convective zone tg^ = dt /dlnTj, 
(dotted line) as a function of the base temperature 7i for pure helium accre- 
tion (X = 0, y = 0.99, Z = 0.01) at ni/ffiEdd = 0.1 and solar abundance accretion 
(Z = 0.71, F = 0.28, Z = 0.01) at m/ffiEdd = 0.1 and 0.01. The convective zone 
reaches a minimum pressure approximately when the equality ^th = %r is first 
satisfied. 



During the burst, convection is very efficient over most of the 
convective region and Nu ^ 1 . Because mixing length theory 
only provides an order of magnitude estimate of Fios.s(3'c), we 
introduce a scaling prefactor A to parameterize the uncertainty 
in its exact value. We thus have Fioss(3'e) = F - Fcomiyc) ~ 
A/^Vad/Vrad = A4flcr/Vad/3K3'c whcrc Tc = T(yc). We as- 
sume A = 1 throughout, though as we show below, jc min is 
sensitive to the value of A. 

In Figure^we show the evolution of the growth time scale 
of the convective zone, fg, = dt/dlnTb, and the thermal time 
scale at the base of the radiative zone, t± = CpTcyc/F\osJjc)- 
During the early stages of a burst, the convective-radiative in- 
terface is located at large column depths and fth ^ fg,. Thus, 
the radiative region cannot thermally adjust to the grow- 
ing convective region and the thermal profile in the radia- 
tive region is unchanged from the pre-ignition profile (see 
[Hanawa & Suaimoto 1982). The initial entropy of the atmo- 
sphere, which is set by Fo = FH + F„a?.u is therefore important 
in determining the convective evolution. ' 

Eventually the convective region reaches a low enough col- 
umn depth that fth = fgr- The radiative flux FiossCjc) can then 
finally diffuse through the entire radiative region without be- 
ing overtaken by the growing convective region. Some of this 
flux will heat the radiative region, while the remainder es- 
capes through the photosphere. The subsequent evolution of 
the convective-radiative interface is determined by the column 

' When the density and temperature first get high enough for helium to 
ignite the burning time scale may be longer than or comparable to the thermal 
time scale. The flux from this early burning may change the pre-burst profile 
of the radiative region slightly, though eventually the burning becomes non- 
linear and the burning time becomes much shorter than the thermal time. 
To examine how this might affect the growth of the convective region we 
computed yc(t) assuming an artificially high pre-burst flux (e.g., 5 X Fq). We 
find that this effect may increase Jcmin by as much as a factor of 3. 



depth at which fth = fgr, i.e., the radiative region continuously 
adjusts to Fioss(3'c), which varies due to changes in the burning 
rate. 

Figure [2 shows results for accretion of pure helium at a 
rate m/zwEdd = 0.1 and solar abundance accretion at a rate 
m/niEdd = 0.01 and 0.1 (corresponding to models HeO.l, 
HHeO.Ol, and HHeO.l as described in Table^. The pure He 
model is similar to the ignition models found by C03 in fits 
to the burst properties of 4U 1820-30. For hydrogen-rich ac- 
creted material, the lower the m the shorter the growth time at 
a given Th- This is because the mass fraction of helium at the 
burning layer, and therefore the triple alpha energy generation 
rate £30, oc Y^, increases as m decreases. 

In models with He-rich burning layers, fgr decreases sharply 
at Ti, > 10'' K due to the sudden rise in energy generation e 
when the '^C(p,7)'^N(Q;,p)'^0 bypass first takes over the re- 
action flow (see When ~ 1 .5 x lO** K^r = fth 10"^ s, 
resulting in millisecond burst rise times (§ 14. 2i . 

In Figure |2l we show the temperature profile at different 
stages for the pure He m/»jEdd = 0.1 model and the mixed 
H/He in/mEdd = 0.01 model. For models accreting hydrogen, 
there is a large compositional contrast between the helium- 
rich matter that is burning and the outer hydrogen-rich ma- 
terial.^ As we demonstrate momentarily, the contrast inhibits 
the outward progress of the convective zone because the burn- 
ing material must get even hotter to become buoyant in the 
overlying hydrogen-rich envelope. Thus, the minimum col- 
umn depth reached by the convective zone yc,min is signif- 
icantly smaller for helium accretors as compared with hy- 
drogen accretors, even when the latter accretes at rh < m^iii 
and thus also bums in a pure helium environment (see also 
Gumming & Bildsten 2000). This difference in convective 
zone evolution between the hydrogen and helium accretors 
influences fth(f) (see Figure ^ and therefore affects the light 
curves during burst rise, as we show in § 14.21 

To appreciate how the compositional contrast influences the 
convective evolution, first consider the ver y ear ly stages of the 
evolution, when fth ^ fgr. As noted in § 12.11 we define the 
convective-radiative interface as the location of neutral buoy- 
ancy, where the density of the radiative solution just equals 
that of the convective solution. The density is therefore con- 
tinuous across the interface, even in the presence of a compo- 
sitional contrast. Hydrostatic equilibrium ensures the pressure 
is also continuous. Thus, to accomodate the jump in mean 
molecular weight across the interface, the temperature must 
increase by a fractional amount {Tc-Tr)/Tc ~ {p^c- IJ-r)/ IJ-c, 
where Tc and jic are the temperature and mean molecular 
weight on the convective-side of the interface, Tr and /x,- are 
those on the radiative-side, and the equality assumes gas pres- 
sure dominates. A thermal wave propagates outward from the 
interface, but only travels a distance « /?(fgr/fth)'''^, where h is 
the scale height at yc, before being outrun by the growing con- 
vective zone. Thus, during the very early stages of the evolu- 
tion, when fth ^ fgr, the time-dependent solution yields a very 
sharp temperature gradient. We approximate the gradient as a 
jump as long as fth > fgr, as shown in Figure |2] Our approxi- 
mation becomes inaccurate as the convective zone extends out 
to lower pressures and fth — > fgr, as then the thermal wave can 
"get ahead" of the growing convective zone. This approxima- 
tion thus prevents us from resolving the very early light curve 

^ Since the particle diffusion time scale is much longer than fgr, there is not 
adequate time for semi-convection to develop. We therefore approximate the 
composition gradient as a discontinuous step function. 
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Fig. 2. — Evolution of the temperature profile for pure helium accretion at rn/rfiEM = 0.1 (left panel) and solar abundance accretion at rh/niEiii = 0.01 (right 
panel) at nine different time.s corresponding to Tj/lO' K = Tq, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5, and 1.7, where Tq is the initial base temperature set by the radiative 
solution with Fq = Fh +f"cmst. The open squares mark the convective-radiative interface Vc. The vertical dashed lines denote the column depth where y = O.Olyi,, 
corresponding approximately to the column depth ywind below which material is ejected by the radiative wind of an RE burst. 



of the burst rise. We resolve the light curve well-before the 
luminosity exceeds the accretion luminosity, however 

Once fth = fgr, there is enough time for heat diffusion (via 
radiation) to smooth out the temperature gradient over a scale 
height, and the radiative region heats up. During this stage 
the convective zone returns to higher pressure and the com- 
positional contrast between burnt and unburnt matter is at a 
fixed location y = jc^min and is thus entirely within the radia- 
tive zone. The discontinuous change in opacity at yc.min results 
in a change of slope of the radiative temperature profile at that 
location, as can be seen in Figure |5] 

We now obtain a rough estimate of jc.min by setting fth = 
fgi and solving for yc given the pre-burst thermal profile set 
by Fi). Using the mixing length theory relation for FiossCjc) 
and assuming an adiabatic profile in the convective zone so 
that Tc = Tb(yc/yb)" yields fth = l>KCpy];'y^-^" /AacXriTl Then, 
since fgr « CpTb/e, we have, upon setting fth = fgr, that 

At yc the density and pressure are continuous and at yc,mm the 
pressure is nearly that of an ideal gas. We can thus elimi- 
nate Tt, from the above expression using Tc = HcTAyc)/ IJ-r and 
T,{yc) ^ OKycFo/acy/'*. Solving for y^ we get. 



(12) 




where the dependence on m is obtained by substituting equa- 
tion QforFzi =Fo/102i ergs cm s . In the lower expres- 
sion, Yt, is the helium mass-fraction in the burning layer and 
we assumed n = 2/5, A = 1, and an energy generation rate 
ei8 = e/lO'^ ergs g"' s"'. For helium burning, cig/Y^ = 5.3 x 
I0^plexp(-4A/Tg)/Ti / Hansen fc KawaleillQQ-l where Tg = 
7/,/10^ K and p(, = ph / 10^ g cm~^. The estimate and its scal- 
ing with burst parameters agrees reasonably well with our full 
evolution calculations presented in § |4] Note the strong de- 
pendence on compositional contrast between the accreted ma- 
terial and the ashes. 

As we show in § |3 the pure helium accretion models and 
the solar abundance accretion models with m < O.OSmEdd 
achieve a super-Eddington luminosity that drives a radiative 
wind capable of ejecting material located at column depths 
y < Jwind — O.Qlyb- Since the convective zone reaches pres- 
sures ycmin ^ Jwind in thcsc modcls, one expects ashes of 
burning to be amongst the wind ejecta. 

3. NUCLEOSYNTHESIS FROM PURE HELIUM BURNING AND 
MIXED HYDROGEN AND HELIUM BURNING 

The full nuclear reaction network used in this work con- 
tains 394 nuclei from H to Sr ranging from the proton drip 
line to a few mass units beyond stability and accounts for 
all relevant proton and a particle induced reactions and the 
corresponding inverse processes. The reaction rates have re- 
cently been updated and a full discussion will be presented 
in a forthcoming paper Here we are only concerned with 
the nucleosynthesis in the early, convective stages of an X- 
ray burst. The relevant nuclear reactions involve only nu- 
clei lighter than iron and we restrict the following discus- 
sion to those nuclei. When no experimental information 
is available, reaction rates are taken from sd-shell model 
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Fig. 3. — Time-integrated net reaction flow for model HeO. 1 {left grid) and model HHeO. 1 {right grid). The thickness of each line indicates the strength of the 
reaction flow relative to the 3a reaction: more than 10% flow {thick solid line), 1%-10% flow {thin solid line), and 0.1%-1% flow {dashed line). Each square 
stands for a proton stable nucleus and filled squares are stable nuclei. Shown are all nuclei below Mn that are included in the reaction network. Note that the net 
flows for (p, 7) reactions can be very inaccurate when both the (p, 7) rate and the inverse (7, p) rate are very large. 
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Fig. 4. — Mass fraction of various isotopes as a function of time {bottom axis) and base temperature {top axis) for model HeO.l. The left panel shows the 
full evolution and the right panel zooms in on the time during which the '^C(p,7)'^N(a,p)'*0 bypass takes over the reaction flow from '^C to '*0. Only those 
isotopes whose peak mass fraction exceeds 0.05 are shown. 



calculations ('Herndl et al."1995X fp-shell model calculations 
(i Fisker et al. 2001), or, the stati stical model NON-SMOKER 
jRauscher & ThielemannI l2000h . Experimental rates were 



taken from the NACRE ("Angulo et al."1999|) and llliadis et"al] 

(2001) compilations with some exceptions. As far as the re- 
action rates relevant to this work are concerned, these ex- 
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whose peak mass fraction exceeds 0.05 are shown. 



ceptions include: the 3 a-reaction where we used the new 
experimental rate from iFvnbo et"ai] (|2005), though for the 
temperatur es relevan t here this rate is identical to the one in 
|Cau2hlan &Fowleij ( 1988). For ^^C(a,"f) we use the reaction 
rate by Buchmann ( 1996). The ^'^0(a,p)^^F reaction is from 
Eahn et al. ( 1996) and the '^N(p,a!)^^0 and ^'*Mg(a,7)^^Si re- 
actions are from Cauahlan & Fowler ( 1988). Weak interac- 
tion rates were taken from . Fuller. Fowler. & Newman (.198 2) 
when a vailable, otherwis e ground state /3-decay rates from 
Nubase JAudi et alJ2003h were employed. Weak interactions, 
however, do not play a role in any of our models. 

We calculate the nucleosynthesis and nuclear energy gener- 
ation only at the base of the convective zone yt,. To estimate 
the integral of energy generation over column depth in the 
calculation of dTb/dt and fgi (see equation |9}, we first take 
the ratio of J e{Th,y)dy to e(Th,yb)yh assuming only a-capture 
reactions. We find that to a good approximation the ratio is 
constant with temperature and we use this constant factor to 
rescale e(Th,yb)yb calculated using the full network. Our ini- 
tial composition is the composition at yi, at ignition and we 
neglect compositional changes due to convective mixing as 
the convection zone moves outward. Such mixing might be 
relevant for the H/He accretion models where the convection 
zone at later stages extends into hydrogen rich regions. How- 
ever, the effect is probably comparable to an increase in ac- 
cretion rate, which also results in a larger hydrogen abun- 
dance. Therefore, our neglect of convective mixing might 
simply translate to a slight underestimate of the accretion rate 
in our sequence of calculated models. 

Figure|3lshows a typical reaction flow during pure He burn- 
ing and Figure |4] shows the composition as a function of 
time and base temperature Tb- After a slow rise from 0.2 to 
0.4 GK temperatures rise quickly, typically within 10-100 ms 
from 0.4 GK to peak temperatures set by yb of 1.8 GK. The 



reaction flow is mainly characterized by a sequence of a- 
captures into the ^^Ar-^Ti region. However, reactions off 
the a-chain also play an important role. For example, be- 
yond Mg, (a,p)-(p,7) sequences carry a substantial part of the 
flow and in fact dominate in some cases, such as for the flow 
from ^^''Ar to '^"Ca, or the flow beyond '^"'Ti. The importance 
of (q;,p)-(p,7) links and their inverses in He burning has been 
jrointed out before (^Hashimoto. H anawa. & SugimotdfT983l 
iTjmmes. Hoffman. & Wooslev 200(F 

The most important effect of off a-chain reactions however 
is the bypass of the '^C(q;,7) reaction, the slowest reaction in 
the a-chain up to ^^S for temperatures > 1 GK, by the much 
faster '^C(p,7)'-'N(a,p)"'0 reaction sequence. The bypass re- 
quires protons as catalysts, which are produced by (a,p) reac- 
tions on heavier a-chain nuclei such as ^"^Mg, ^^S and ^''Ar. 
Initially, He burning proceeds via the 3a reaction, building up 
'^C with some weak reaction flow through '^C(a,7) and suc- 
cessive a captures. At ^"^Mg, ■'^S, and ■'^Ar, (a,p) branchings 
are significant and lead to the build up of a small amount of 
protons. For example, in the pure He m/iriEdd = 0.1 model the 
proton abundance reaches 3x10"^ at a temperature of about 
1 GK, which is sufficient for the '2C(p,7)'^N(a,p)"'0 bypass 
to take over the reaction flow from C to '^O. At this point 
a runaway effect sets in as any increase in flow beyond '^C 
increases the production of protons, which further increases 
the efficiency of the '^C(p,7)'^N(a,p)'*0 bypass. Only 2 ms 
later at a temperature of 1.2 GK the flow from '^C to '^O ex- 
ceeds the flow via the 3a reaction and '^C is quickly depleted. 
As a result, the '^C(p,7)'-'N(a,p)'^0 reaction sequence com- 
pletely dominates the reaction flow towards heavy elements 
during the convective burning stage. This is a very robust ef- 
fect, which is independent of the initial metallicity, and will 
always occur during He burning under high temperature con- 
ditions that lead to the production of^^Mg or heavier nuclei. 

The rapid onset of the '^C(p,7)'''N(a,p)'^0 bypass leads to 
a sudden spike in energy production, which causes the drop 
in fgr and t± to ^ ms time scales and the second dip in Vc, 
as seen in Figures [J and |2j respectively. As we show in §|3 
the bypass occurs in all models burning helium-rich material 
regardless of m (see Figure |6jl, and is directly observable as 
a rapid increase in the rise time of the burst light curve (see 
Figure |8}. For comparison, we also calculated the pure He 
m/riiEdi = 0.1 model with a simple a-chain network, the re- 
sults of which are shown in Figure |6l |8] and Another 
consequence of bypassing the '^C(a,7) reaction is the deple- 
tion of '^C and the enhancement of the production of heav- 
ier a-chain nuclei in the final composition. As Figure [T3l 
shows, a pure He nt/rh^ii = 0.1 model with just an a-chain 
network would predict the main product to be '^C, while with 
the full reaction network the main products are ^^Si and ^^S 
with negligible amounts of carbon. In addition, the presence 
of protons leads to proton capture reactions producing small 
amounts of nuclei off the a-chain. The very small amounts of 
nuclei around A = 60 are produced by proton capture of iron 
peak nuclei present in the accreted material. 

For pure He accretors, except for model HeO.lb (see Ta- 
ble 0, the composition at the time the convective zone has 
retreated back to the base is mainly doubly magic '*'^Ca and 
unburned "^He, with some ^^Ar and "^"^Ti. The most abundant 
non a-chain nucleus is ^^K. HeO. lb ignites deeper, at a higher 
mass density, resulting in a higher peak temperature. Burning 
therefore proceeds beyond '^'^Ti with a sequence of alternating 
(a,p) and (p,7) reactions up to ^^Fe. Besides unburned He, the 
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composition is then dominated by '*'*Ti and "^^Mn. For mixed 
H/He accretors that undergo pure He burning at low accretion 
rates (e.g., models HHeO.Ol and HHeO.OlX) the composition 
is very similar to the pure He accretion models except for a 
somewhat lower '"'Ar and higher ''^Ti production. This dif- 
ference is due to the faster rise in temperature in the H/He 
models, a consequence of the smaller convective extent in the 
presence of hydrogen and thus a smaller convective mass to 
heat up for the same energy input (see fgr curve for HeO. 1 and 
HHeO.Ol in Figure[lli. 

Though some protons are produced during He burning, the 
presence of large amounts of hydrogen at ignition in models 
such as HHeO.l {X = 0.4) and HHeO.OlXZ {X = 0.08) leads 
to a drastic change in the nuclear reaction sequence. Figure|3] 
shows a typical reaction flow and Figure |5] shows the com- 
position as a function of time and base temperature Tb- The 
flow is governed by the ap-proces s, a sequence of (a,p) and 
(p,7) reactions beginning at ( Wallace & Woosle\ill981l) . 
For the larger proton abundances in HHeO. 1 several parallel 
side branches involving more neutron deficient nuclei are es- 
tablished. The main product at the time the convective zone 
returns to the base tends to be ''**Ca (37% by mass) together 
with small amounts of "^^Ti (7% by mass). Since He is abun- 
dant during the convective burning phase and since the burn- 
ing is largely confined to Z < 40 nuclei, (a,p) reactions are al- 
ways faster than /3-decays. Thus, the rapid proton capture pro- 
cess (rp-process) of Wallace & Wooslev (1981), a sequence 
of proton captures and /3-decays, does not yet play a role at 
this early burst stage. It will take over at later times, pro- 
ducing post-peak light curves like that seen in GS 1826-24 
iOallowav et al..2004.) . 

4. RESULTS 

The radiative wind generated by the super-Eddington lu- 
minosity of an RE burst will eject some ashes of burn- 
ing if: (1) the fraction of the accreted mass ejected during 
the wind satisfies AMw/Macc = ywindM > yc,mm/yb, where 
Mace AirR^yt « 4 x lO^' g (y^/S x 10** g cm'^), and (2) the 
wind is generated subsequent to the convective region reach- 
ing jc < jwind- In § 14. II we calculate yc^min and demonstrate 
that condition (1 ) is satisfied for a wide-range of burst igni- 
tion models. In § l4.2l we calculate the light curve during burst 
rise and show that condition (2) is satisfied for those same 
models. 

We first note that since the gravitational binding energy 
at the surface is « 100 times greater than the helium burn- 
ing energy release per unit mass, at most w 1 % of the atmo- 
sphere can be ejected by the wind. To obtain an estimate of 
AMw/Macc, suppose the surface luminosity as measured by 
an observer at infinity L is super-Eddington so that the energy 
loss rate at the surface E ~ L + M^.c-fjl+ z)- At the photo- 
spheric radius of the wind rph ^ R (see ^5.3> the luminos- 
ity is very nearly the local Eddington value Lph ~ iEdd('"ph) = 
47rGMc(l +Zph) /k ~ AttGMc/k. Equating the energy loss rate 
at the surface to that at the photosphere yields the mass-loss 
rate due to the wind 
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4.1. Dependence of the convective extent on burst parameters 

In Figure |6l we show the time variation of the convective 
extent jc during a burst, with the progression marked by the 
temperature at the base of the burning layer, 7i, rather than 
by time. We consider burst ignition models assuming both 
mixed hydrogen/helium accretion and pure helium accretion 
(see Table[0; the latter models are similar to those obtained by 
C03 in fits to the burst properties of 4U 1820-30. We consider 
models over a range of accretion rate, ignition column depth, 
metallicity, and crustal heat flux. 

The evolution of jc is sensitive to the assumed ignition con- 
ditions and to the energy production from non-a-capture re- 
actions (see §|3}. Ignition conditions that maximize the peak 
flux of a burst and minimize the entropy of the atmosphere 
during the fuel accumulation stage yield bursts with the small- 
est yc, min- This is because the larger the peak flux, the faster 
Tt, rises, and hence the shorter the growth time scale fgi at a 
given Th- The convective region must therefore reach out to 
lower column depths before fth(3'c) = fgr- As illustrated in Fig- 
ureQ such ignition conditions are best satisfied at low m. The 
reason is twofold. First, for rh < 0.04mEdd and solar abun- 
dance accreted material a pure helium layer develops at the 
base of the atmosphere and the peak luminosity is consider- 
ably greater than that of bursts with mixed hydrogen/helium 
burning layers. Second, because the pre-burst flux oc m, the 
lower rh the smaller the initial atmospheric entropy. 

For the models considered, the fraction of accreted mass 
lying above the convective region at its maximum extent is in 
the range 10""* < AM /Mace < 10"^ making it likely that ashes 
get ejected. 

4.2. Light curve during burst rise and radiative winds 

We determine the rising portion of the burst light curve by 
calculating the radiative flux loss at the photosphere Fph as 
Ti, increases. To obtain Fph we integrate equation ([3 inwards 
(assuming the radiative zero solution at the outer boundary), 
varying the flux at the top until the radiative solution intersects 
the convective solution at the column depth jc- Initially, fth > 
fgi- and Fph = 0. Eventually fth = fgr and Fph rises as the radiative 
region heats up due to the rise in Fioss(3'f)- In Figure Owe 
show the rising portion of the light curve for the same models 
as Figure Hand also for a pure helium model using only the 
a-chain network. We plot the luminosity L = AirR^F-pk as a 
function of f - fo, where fo corresponds to the time when the 
equality fth = fgr is first satisfied. 

The rise time at a given time f is set by the thermal time at 
the convective-radiative interface fth(3'e[f])- As the differences 
between the HeO.l and HeO. la models illustrate, the shape 
and rise time of the light curves are sensitive to the energy 
generation rate. The sharp sub-ms rise in models HHeO.Ol 
and HeO.l occurs at the onset of the '^C(p,7)'-'N(Q;,p)'*0 



9 



CM 
I 



E 
o 



> 
c 
o 
u 



10^ h 



10^ k 



10^ hr 



10^ 



10^ 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




HeO.OI 



HeO.lb 



HeO.2 = HeO.IQ 

_l I I I I I I I I I I I I I I I I I I I L_ 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

HHeO.18 

HHeO.08/ 




HHeO.01 3 
HHeO.OIX 



:: HHeO.OIXZ 



J I I I I I I I I I I I I I I I I I I I L_ 



0.5 



1.5 2 0.5 

Base Temperature (10^ K) 



1.5 



Fig. 6. — Evolution of the top of the convective zone Vr as a function of base temperature T/, for burst models with pure helium accreted material (left panel) and 
mixed hydrogen-helium accreted material {right panel), as described in TablefTI The two unlabeled helium models are, from bottom to top, He0.04 and HeO. 1. 
Model HeO. la {dash-dot line) was calculated using a simple a-chain reaction network. The horizontal dashed line denotes the column depth where v = O.OlVj,, 
con'esponding approximately to the column depth y^ind- 
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Fig. 7. — Minimum column depth reached by the convective zone yc^mm as 
a function of m in units of the Eddington rate. Results are shown for the pure 
hehum accretion models He[ffi] and the mixed hydrogen/helium accretion 
models HHe[m]. The horizontal dashed line denotes the column depth where 
y = 0.01>'j,, corresponding approximately to the column depth y„ind- 



bypass (§ 1^. The bypass leads to a spike in energy pro- 
duction which in turn drives the convective zone to lower 
pressures. Since T{yc) = is very high by this time, fth = 
CpTcyc/Fjossiyc) = 'inCpy^/AacXnT^ becomes very small (see 



Fig. 8. — Rising portion of the burst light curve for the same three mod- 
els shown in Figure0(,so/iy lines). The dash-dot line is the light curve for 
the HeO. la model calculated using just an a-chain reaction network rather 
than the full network. The luminosity is plotted as a function of the offset 
time t — tfi where /q corresponds to the time when the radiative region first 
begins to evolve away from its pre-burst state (i.e., when the equahty fth = fgr 
is first satisfied). The dotted lines denote the constant accretion luminosity 
for the models and the dashed lines denote the Eddington luminosity at the 
photosphere. 



Figure n). In the HeO.l model, the light curve initially rises 
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Fig. 9. — The local Eddington luminosity Lgdci as a function of column 
depth at the end of the calculation, when the convective zone has fully receded 
back to the base. The discontinuity at Vcmin for the H/He models is due to a 
discontinuity in hydrogen abundance between burnt and unburnt material. 



Fig. 10. — Rise time of the burst light curve as a function of rh in units of 
the Eddington rate for the same models as FigurelTl Two sets of rise times are 
computed for each model: rdctcct (solid lines) is the instantaneous rise time 
dt / d In L when the luminosity first exceeds the accretion luminosity, corre- 
sponding to when an observer can first detect the burst; T|^2 (dash-dot lines) 
is the time it takes the burst to rise from 10% to 50% of its peak luminosity. 



relatively slowly (> 10 ms) and then, once the temperature is 
high enough for the bypass to take over, it sharply rises. This 
is distinct from the HHeO.Ol model which rises sharply at the 
outset because f = fo nearly coincides with the onset of the 
bypass. 

Depending on the ignition conditions, the base tempera- 
ture can get sufficiently high that the flux becomes super- 
Eddington. In Figure |9l we show the local Eddington lumi- 
nosity LEddiy) as a function of column depth for the mod- 
els HHeO.Ol, HHeO.l, and HeO.l at the end of the calcu- 
lation, when the atmosphere is completely radiative (i.e., 
yc = yb)- The increase in LEdd(3') with depth is the result 
of the decrease in Kes with temperature due to relativis- 
tic corrections. The discontinuity in LEddiy) for the mixed 
H/He models occurs at yc.min and is due to the discontinu- 
ity in hydrogen abundance between the burnt and unburnt 
material. Such a discontinuity could explain the possible 
bimodal distribution of the ob served burst peak luminosi- 
ties (see iKuulkers et al.l '2003 for a recent discussion), as 
first pointed out by [Sugimoto, Ebisuzaki, & Hanawa ( 1984^ 
and investigated in greater detail by Ebisuzaki & Nakamura 
0988). They proposed that the ejection of the hydrogen-rich 
envelope during radius expansion burst occasionally exposes 
the underlying helium-rich region and results in the Edding- 
ton limit transitioning from that of hydrogen-rich material to 
that of hydrogen-poor material. 

We find that the peak luminosity exceeds LEddiy) for m < 
0. ImEdd if the accreted material has solar abundance {m can be 
higher if the material is helium-rich) and that the Eddington 
value is first exceeded near the top of the atmosphere. Such 
systems are expected to develop a radiative wind. Eventually, 
the flux loss exceeds the flux from nuclear burning and the at- 
mosphere cools, though we do not calculate this portion of the 
light curve. Note that since the convective zone reaches a min- 
imum pressure at f ~ fo, bursts that become super-Eddington 



achieve yc < jw ind well before the wind is generated. 

In Figure^] we show the dependence of the rise time on 
m for the pure He models and the mixed H/He models. We 
calculate rise times during the early portion of the rise, when 
the surface flux just exceeds the accretion flux (raetect), and 
during the later portion of the rise, taken to be the time for the 
burst to rise from 10% to 50% of its peak luminosity (ti/2)- 
Since fth(3'e.min) sets the initial rise time, models with smaller 
values of ycmin/^c^min ^^^'^ ^avc Smaller Tdetect. which, to a 
first approximation, corresponds to models with lower m. For 
helium accretors, ti/2 is nearly independent of m because the 
convective zone evolution at late times is very similar among 
models that differ only in m (see Figure |6}. By contrast, for 
accretion of hydrogen-rich material, Ti/2 increases by a fac- 
tor of w 50 between m/niEdd = 0.02 and 0.1. In their mul- 
tizone numerical simulations of X-ray bursts, 'Woosle v et alJ 
(2004) consider models that accrete mixed hydrogen/helium 
with solar abundance at m/mEdd = 0.02 and 0.1. They fol- 
low several sequences of bursts for each model and obtain 
values of T1/2 in the range 0.51 x 10~^-32.1 x 10~^ s for 
ni/riiEdd = 0.02 and 0.51 -0.66 s for m/riiEdd = 0.1, both in 
good agreement with our estimates. Our rise times are also 
in broad agreement with t he relevant models of earlier nu- 
merical simulations (e.g.. lTaamlll98U lAvasli & JossI 119821 
■Wallace . Wooslev, & WeaveJ 19821) ! 

5. DETECTING THE NUCLEAR BURNING ASHES 

In § |4] we showed that the convective zone reaches suffi- 
ciently low pressures during RE bursts that yc^min ^ Jwind over 
a broad range of burst parameters. Thus, ashes of nuclear 
burning can be ejected by the radiative wind of RE bursts. 
Furthermore, when the photosphere settles back down to the 
NS surface following the RE phase, it is laced with heavy- 
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element ashes. In this section we desc ribe the composition 
of the ejected and exposed ashe s (§ 15. 1> and address whether 
they can be detected. In § l5.2l we discuss the possibility that 
some of the p-nuclei found in the solar syste m ow e their ori- 
gin to ash ejection during RE bursts. In §§ 15. 31 and I5.4l we 
determine the expected strength of photoionization edge fea- 
tures from ashes ejected in the wind and those exposed at the 
NS surface. 

5.1. Composition of ejected and exposed ashes 

Just prior to the onset of the wind, the convective zone has 
receded to the base and the atmosphere has a stratified compo- 
sitional structure. Throughout the region y < yc.min, the com- 
position is that of the un proc essed accreted material. As we 
show in Figures ^3 and 1121 for y > yc^mm the composition 
is determined by the burning stage at the moment ydt) = y 
during the convective zone's retreat to the base. Results are 
shown for models HeO. 1 and HHeO. 1 . The discontinuity at 
y = 3 X 10^ g cm"^ in model HeO.l is due to the rapid nu- 
cleosynthesis that occurs once the '^C(p,7)'-'N(Q!,p)'^0 by- 
pass takes over the reaction flow from '^C to '^O (§|3}. The 
composition at the base of the wind is approximately that at 
y = O.Olyb = 3 X 10^ g cm"^. In Figure [O] we show the ash 
composition at y/yh = 0.01 ~ AM^/Macc for models HeO.l, 
HeO. la, and HHeO.l. While helium comprises the largest 
fraction of the mass in all three models, the models' overall 
abundance distributions differ significantly from one another 
due to the impact of the '^C(p,7)'-'N(a,p)'^0 bypass and dif- 
ferences in the initial composition at the base of the burning 
layer. Some of the proton-rich isotopes shown in the above 
figures have half-lifes comparable to the duration of the atmo- 
sphere's convective phase. Though we do not show the decay 
products in these figures, we do account for such decays in 
our calculation of photoionization edges below. 

5.2. Ejection of p-nuclei 

In their n umeric al simulations of X-ray bursts, 
IWooslev et al.l i2004ft find that bursts ignite in the ashes 
of previous bursts. If so, these endpoint ashes, which are 
heavier than those processed during the ongoing burst at 
the time yc/yi, = 0.01, are mixed throughout the convective 
region (we do not show such ashes in Figures [T2l orll3>. 
Using a reaction network that extends up to Xe, Schatz et al.l 
(jlOOl) find that in models where solar abundance material 
is accreted at m = 0.3;7tEdd (i-C-, burning in a hydrogen-rich 
environment), the endpoint of rp process burning is a closed 
SnSbTe cycle that naturally limits rp process nucleosynthesis 
to light p nuclei. They find overproduction factors (relative to 
solar abundance) of ~ 10*^ for the p nucle i ^^Mo and ^^Ru and 
- 10'' for the p nucleus ''**Ru (see also ISchatz et alJll998h . 
Standard p-process scenarios are unable to adequately explain 
the observed solar system abundances of these p-nuclei (for a 
review, see Wallerstein et al.^1992) . 

Whether ash ejection during RE bursts can account for the 
observed solar system abundances depends on the amount of 
p-process material ejected into the interstellar medium by an 
RE burst and the event rate of such bursts over the Galaxy 
lifetime. In order to produce large amounts of p-nuclei, the 
burning layer must be hydrogen-rich (m Ki 0.3mEdd) while 
RE bursts require helium-rich burning layers and thus m < 
0.05mEdd- However, this does not preclude RE bursts from 
ejecting p-nuclei, as accretion rates in bursting low-mass X- 
ray binaries are observed to vary by factors of a few, with in- 
dividual systems undergoing transitions from hydrogen-rich 



to helium-rich burning over year time scales (as evidenced 
by variations in burst duration and peak fluxes; see e.g., 
ICornelisse et alJl2003h . Another limitation to how much p- 
nuclei an RE burst can eject is that the p-nuclei tend to be 
produced at the bottom of the burning zone where the hydro- 
gen is consumed. Thus the ignition and base of the burning 
zone may lie well above the p-rich ashes of a preceding burst. 
Unless the convective zone dredges up ashes from deep be- 
low, the amount of p-nuclei mixed into the convective zone 
may be small. 

To determine the fractional amount rj of ashes that must 
be ejected in order to account for the observed solar sys- 
tem abundance of p-nuclei, assume a p-nuclei overproduc- 
tion factor £ = 10^ and a galax y disk mass Mjisk = 4 x lO"' 
('Klvp in. Zhao. & Somervillel2002l) . Currently, there are « 10 
active X-ray burst systems at M ~ 10~^Mq yr"' and « 100 at 
M - IO-'^Mq yr-i (iLewin. van Paradiis. & Taamlll995l) . If 
we assume this is representative of the population count over 
the galaxy lifetime, then the total amount of mass accreted by 
all RE burst systems over 10'° yr is Macc.tot ^ lOOM©. A frac- 
tion ?7 = 0.4(Mdisk/4 X 10"'Mo)(e/109)-i(Ma,c,tot/100MQ)-i 
of all accreted material must therefore be ejected. Thus, 77 is a 
factor of ~ 100 too high given that only « 1% of all accreted 
matter is ejected. The discrepancy can be overcome if, for ex- 
ample, the Galactic distribution of p-nuclei is inhomogeneous 
(i.e., the solar abundance is higher than the Galactic mean by 
a factor of ^ 100) or the population of bursting systems was 
much larger at earlier times in our galaxy. 

5.3. Spectral Edges in Wind Outflow 

During the RE phase, an optically thick, transonic, 
radiation-driven wind forms. The sonic point of the wind 
lies 10- 100 km above the NS surface and the photo- 
sphere, defined as the location where the effective opti- 
cal dept h Tj, = Kpr is near unity, is a factor of ^ 10 far- 
ther out ("Paczvuski & Pr oszvnskilll986t Lloss & Meiialll987t 
Nobili. Turolla. & Lapidusl ll994r Once matter reaches the 
sonic point it is essentially unbound from the NS, and is 
ejected to infinity. 

The ejected matter is at sufficiently low temperature during 
the RE phase that some heavy elements bind with one or more 
electrons. The resulting column density of hydrogen-like 
ions above the photosphere is thus A^wind ^ f(A,Z)C,{T, p)Ne — 
10^° cm"^(C//10~'*), where /(A,Z) is the abundance by num- 
ber of element Z with mass number A, C{T,p) is the fraction 
in the hydrogen-like state at a given temperature and den- 
sity from Saha equilibrium, and Ng ~ cr^J is the electron col- 
umn density. This corresponds to an optical depth to a pho- 
toionization edge of T w A^windfbf ~ 1, where the bound-free 
cross section crbf(-£') = 6.32 x 10 "'^ cm^(Ee/E)^Z~^ and Eg c f. 
13.6Z^eV is the edge energy (iRvbic ki & LightmanllT979l) . 
The effective temperature is low enough at the photosphere 
that even for metals at solar abundances r ~ 1. That the 
ejected ashes have abundances much larger than solar im- 
proves the likelihood of detection. 

We now calculate A^wind more exactly and determine the re- 
sulting equivalent width (EW) of the photoionization edge. 
The column density of hydrogen-like ions for a given species 
is 

A^wi„d = ^^^/ aT(r),p(r)]p(r)dr. (15) 

The recombinations are nearly instantaneous on the time 
scale of the wind flow, i.e., fi-ec ^ fflow ~ r/v{r) « 0.1 s. 
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Fig. 11. — Ash composition as a function of column depth after the convective zone has completely receded (i.e., yc =yb)- Results are shown for models 
HeO. 1 {left panel) and HHeO.l (right panel). Only those isotopes whose peak mass fraction exceeds 0.01 somewhere in the range yr.min < .V < are shown. The 
vertical dashed line denotes y = 0.0 lyj, ~ jwind- 
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Fig. 12. — Same as Figure lTTI exceDt only those isotopes with A > 50 whose pealc mass fraction exceeds 5 X 10 somewhere in the range yc.mm ^ y ^ yb ^e 
shown. 



where v(r) is the fluid velocity at r and Uez = l/nc(CTfbv) « 
QmZ-^p-_lT~^'^B{ZJ)&, with p_6 = /9/10-«gcm-^ 
T^ = T/\QP K, B(ZJ) = expW /^°°exp(-f)c/lnf, x = EjkT, 
and (Tfh related to through the Milne relation 

jRvbicki & Lightman 1979). We use the results of 
IPaczvn ski & Proszvnski (1986), who calculate general 
relativistic models of radiation-driven winds from NSs, 



to obtain values for rph, r(rph), and jo(rph) as a func- 
tion of the mass outflow rate Mw For example, at 
Mw = 10"^ g s-i, rph = 3 X 10^ cm, Tph = 5 x 10^ K, and 
/Oph = 5 X 10"^ g cm"^. The wind duration is longer 
than fflow so we a ssume that Mw — A-kt^pv (s e e e.g. , 
I.T0SS & Melial 11981) . IPaczvnski & Proszvnskil (Il986h 
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Fig. 13. — Composition of material processed during a burst at the time 
the convective zone is located at jc = O.OIy;, ~ Vwind 2nd is receding to higher 
pressures. Results are shown for models HHeO.l, HeO.l, and HeO.la. 




Fig. 14. — Equivalent width of the photoionization edge due to ashes 
ejected in a wind as a function of the mass outflow rate for the models HeO. 1 
{solid lines) and HHeO.l {dashed lines). 



show that the velocity is nearly constant beyond the 
photosphere, i.e., v(r > r^\) Vph ^ 10^ cms"', so that 

The gas above the photosphere is Compton-heated by the 
hot photons originating in the photosphere. Since the Comp- 
ton heating time scale ^ kTnieC- /E^Fajh ~ 10"^ s ^ tsow, 
where is the photon energy and F the flux, the gas temper- 
ature is nearly constant out to radii well above the photosphere 
(poss & Melia 1987; Nobili. Turolla. & Lapidus 1994). For 

^ '"ph, tc > fflow and the gas cools adiabatically, though 
the density in this region is so low that in calculating A^wind 
we assume T{r > rph) = Tph. For the more abundant species 
shown in Figure we obtain column densities in the range 
A^wind--10"'-102i cm-2. 

To determine the EW of the photoionization edge we inte- 
grate over the optical depth above the edge, = A^wind<7bf(£). 
Thus, assuming an effect ively cold atmosphere, EW^ = f l 1 - 
exp[-N^M<ybf(E)]}dE (B ildsten. Chang. & Paerelsl2003h . In 
Figure [21 we show EW^. as a function of for some of the 
ejected ashes of models HeO.l and HHeO.l. In calculating 
the EWs we accounted for the decay of those isotopes whose 
half-lifes are shorter than 1 s, and therefore shorter than the 
wind duration. Thus, 24si^24j^g^ 285^285^^ ^'^Fe^'^'^Ci, and 
62(3e^622n^ ah from model HHeO.l. As increases, Tph 
decreases and the hydrogen-like fraction increases. For the 
heavier ashes (e.g., *°Zn and ^^Zn) I even at low M„ so 
that A^wind(-Z ^ 30) and EWe(Z ~ 30) are essentially set by 
f{A,Z), as noted above. The computed values of EW^ are 
within the range accessible by current X-ray telescopes. 

Several observational studies report the possible pres- 
ence of spectral lines and edges during the RE phase 
of bursts. Indications of such featu res were found in 
RXTE/PCA observations of GX 11+2 (Kuulkers et al. 2002), 
4U 1 722-30, 4U 2129H-12, XB 1745-25 iKuulkers et al.i 
l2(K)l . and 4U 1820-30 (iStrohmaver & Brownl 120021 



IKuulkers et alJ 120031), in Ginga/ LAC observations of 4U 
2129H-12 (wan Paradiis et al.'1990'), in EXOSAT observations 
of EXO 1747-21 (Magnieretal. 1989), and Tenma obser- 
vations of 4U 1636-53 (Waki et al. 1984). In general, these 
studies report significant residuals in black-body spectral fits 
that cannot be explained by spectral hardening alone. Higher 
resolution observations with Chandra and XMM-Newton are 
clearly needed. 

5.4. Spectral Edges from NS surface 

Following the burst peak, the atmosphere cools, the flux 
becomes sub-Eddington, and the wind turns off. The pho- 
tosphere, still laced with heavy-element ashes, settles back 
down to the NS surface. As in the case of spectral edges 
in the wind outflow, whether these surface ashes can be de- 
tected depends on the column density, N^mi, of ashes that 
are not fully ionized. Here A^surf — fiA,Z)C,{T , p)y^\i/ nip, 
where jph « 1 g cm"^. As the surface temperature T de- 
creases, A^surf increases, although edges cannot be detected 
once L = AirR^aT'^ < Lace, where Lace is the accretion lu- 
minosity. The decay time scale of RE bursts (i.e., the 
time during which Lace < L < Ledd) is typically around ^ 
10 s. The downward drift speed of a nucleus with Z ~ 30 

in a pure H atmosphere is v w 1 cm s~^Tj^^ / p so that it 

takes « 1 s()'/l g cm~2)/^^^^^ for such a nucleus to fall 
a column depth v ("Bildsten. Salo eter. & WassermanI Il992t 
Bildsten. Chang. & Paerels 2003). Thus, the residence time 
of ashes in the photosphere also limits the detectability of 
edges, as will high NS rotation rates (.Ozel & Psaltis 20031 
lOiang. Bildsten. & Wasserman 200^. 

In Figure^lwe plot EWg = J{ 1 - exp[-A^surfO'bf (L)] } dE as 
a function of L/ Lace for the same models as Figure[2] For ref- 
erence, LEdd/i'ace — 10 and 30 for models HeO.l and HHeO.l, 
respectively. The values of EW^ are again within the range 
accessible by current X-ray telescopes and a measurement of 
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Fig . 15. — Equivalent width of the photoionization edge due to ashes resid- 
ing in the photosphere at the NS surface for models HeO. 1 (solid lines) and 
HHeO. 1 (dashed lines). The EW is plotted as a function of the ratio of surface 
luminosity to accretion luminosity for a cooling NS atmosphere following the 
RE phase. Detecting such edges might allow for a measurement of the NS 
gravitational redshift. 



the NS gravitational redshift may be possible. 

6. SUMMARY AND CONCLUSIONS 

We have shown that during a radius expansion Type 1 X-ray 
burst the ashes of nuclear burning can be ejected by the burst's 
radiative wind. Specifically, we solved for the evolution of 
the atmosphere's thermal structure and found that in sys- 
tems accreting pure helium, such as 4U 1820-30, and in sys- 
tems accreting mixed hydrogen and helium at m < 0.05;7jEdd, 
the convective zone reaches sufficiently low pressures that it 
lies within the wind region located at pressures < 1% that 
at the base. Hence ashes of burning, mixed throughout the 
convective zone, can be amongst the wind ejecta. Previ- 
ous studies have also found that the convective zone can 
reach < 1% the base pressure for low m bursts (Jossl|1978[ 
[Taam 1981; AvasH & Joss 1982; Hanawa & Sugimotoi |T982 ; 
iHanawa & Fuiimoto 1984; Wooslev et al. 2004). However, 
they focused on numerical simulations of bursts and typically 
did not resolve the low-pressure zones convection reaches. 
Furthermore, such studies only explored a limited range of 
burst parameter space. Our analytic treatment enabled a sur- 
vey of the dependence of the minimum pressure reached by 
the convective zone on a variety of burst parameters such as 
m, the composition of accreted matter, and the pre-burst ther- 
mal state of the atmosphere. We have compared some of the 
results of our analysis, such as the burst rise times, to those of 
numerical simulations and found good agreement. 

To calculate the nuclear energy generation rate and ash 
composition we used an updated full reaction network that 
contains 394 nuclei and all relevant proton and a particle 
induced reactions and the corresponding inverse processes. 
We have found that in models where burning initiates in a 
pure helium layer (i.e., those systems accreting mixed hydro- 
gen and helium at m < 0.04mEdd or systems accreting pure 



helium), protons are produced by (a,p) reactions on heav- 
ier a-chain nuclei such as ^'^Mg, •'^S, and ^*Ar. The proton 
abundances achieved are high enough to enable the bypass 
of the relatively slow '^C(a,7) reaction by the much faster 
'^C(/?,7)'-'N(a,/9)"'0 reaction sequence. The bypass leads to 
a burst in energy production that pushes the convective zone 
to low pressures and results in a rapid, ~ 10"^ s, rise in the 
burst light curve (see Figure [Tot. The bypass also enhances 
the production of heavier a-chain nuclei, resulting in heavier 
ejected ashes. 

For specific burst models we determined the composition 
of ejected ashes and calculated the expected column density 
of hydrogen-like nuclei using models of relativistic radiation- 
driven winds. We then computed the EW of the photoioniza- 
tion edge for the more abundant hydrogen-like nuclei in the 
wind. We carried out a similar procedure to determine the EW 
for those ashes that remain bound to the NS and thus reside 
in the photosphere after it settles back down to the NS sur- 
face. We found EWs in the range 10- 1000 eV (see FiguresIT^ 
and I15> suggesting that the edges can be detected with cur- 
rent high resolution X-ray telescopes. Detecting them would 
directly probe the nuclear burning. Those edges formed at 
the surface may also provide a measurement of the NS grav- 
itational redshift and thus help constrain the NS equation of 
state. 

If bursts ignite in the ashes of previous bursts then some 
of these processed ashes, which are thought to contain large 
overabundances of p-nuclei relative to solar, are also ejected 
in the wind. We showed that at least ^ 1% of the p-nuclei 
observed in the solar system may originate in X-ray bursts. 

We did not account for the affect of a laterally spreading 
burning front during burst rise nor the influence of rotation, 
though these may significantly alter the convective structure 
and therefore the conditions under which ashes are ejected. 
Bursting systems have magnetic fields B < 10'' G, so the mag- 
netic energy density, B'/^tt, is considerably smaller than even 
the minimum convective energy density pVconv — gycmin — 
5 X 10'** erg cm""*. Magnetic fields are therefore unlikely to 
significantly affect the convective evolution. 

Although low m systems accreting a mix of hydrogen and 
helium yield X-ray bursts with convective zones that reach 
low pressures, their burst recurrence times are long and of- 
ten irregular making them difficult targets to monitor given 
a narrow window of observing time. More promising are 
systems in which the neutron star accretes helium-rich ma- 
terial (Xiis « 0.9- 1) from an evolved companion, such as a 
cold helium white dwarf. The binary 4U 1820-30 is thought 
to reside in such a system as evidenced by its ultracompact 
nature (Porb = 11.4min; Stella, Priedhorskv, & White 198^) 
and the fast rise times, decay times, and a-values observed 
during bursts (Gumming 2003, and references therein). The 
system also exhibits radius expansion bursts with fairly reg- 
ular burst recurrence times of on ly a few hours (see e.g., 
GorneHsse et al. 2003, Kuulkers et alJ2003h . 

Observations of four candidate ultracompact binaries have 
shown an unusual Ne/O abundance ratio in the absorption 
along the line of sight, w ith ratios several times the interstel- 
lar medium (ISM) value ( Uuett. Psaltis. & Ghakrabartvll200U 
Juett & Ghakrabartv 2003). Two of the systems, 4U 1543-624 
and 4U 1850-087, have shown variations in the Ne/O ratio 
in follow-up observations (Juett & Ghakrabartv 2003, 2005), 
suggesting either variations in the ionization state of the Ne 
and O or variations in the intrinsic abundances. Both pos- 



15 



sibilities imply the absorption is due to material local to the 
binaries. 

One explanation for the unusual ratios is that the intrin- 
sic abundance of Ne and O is the same as the ISM but the 
O is in a higher ionization state than Ne, leading to an ap- 
parent enhancement of the Ne/O ratio. Another possibility 
is that the degenerate donors in these ultracompact binaries 
have Ne/O abundances above the solar value. In particular, the 
donors may be the chemically fractionated cores of C-O-Ne 
or O-N e-Mg white dwarfs that have previously crystallized 
dSchulz et al. 1001; Juett & Chakrabartv 2003; Sid oH et all 
2005ir or, in some cas es, the central mass of a He white dwarf. 
in 't Zand et alJ (I2005h have shown that the latter donor type 
is consistent with both the properties of the long X-ray burst 
observed from the ultracompact binary 2S 0918-549 and that 
systems enhanced (and constant) Ne/O ratio. 

A third possibility we now propose is that the donors in 
these ultracompact systems are He white dwarfs, as in 4U 
1820-30, and that the accretion of helium rich matter results 
in RE bursts and ejection of ashes with highly non-solar abun- 
dances of Ne and O. Radius expansion bursts have indeed 
been seen in two bursts fr om 2S 0918-549 (Cornelisse et al. 
I2002tlin''t Zand et alJ2005ft . Though only a single realization, 
the Ne/O ratio by number at y = jwind for model HeO. 1 shown 
in FigureE]is ^ 1 .0, comparable to that seen in the four sys- 



tems. In order for this explanation to work, the ejected ashes 
must either somehow remain in the environment of the binary 
or be continuously replenished by periodic RE bursts. The 
variations in the observed Ne/O ratio in observations of 4U 
1543-624 and 4U 1850-087 could just reflect the time since 
the last RE burst. 

Finally, we note that ash ejection may also occur during 
superbursts that undergo photospheric radius expansion. The 
amount of mass ejected by superbursts may be much larger 
than ordinary RE bursts, suggesting even higher column den- 
sities of ejected ashes. Furthermore, the ejected ashes are 
likely to be heavier since the carbon that fuels the superburst 
lies underneath the heavy rp-ashes of normal bursts. 
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